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Abstract. Previous models of the combined growth and migration of protoplanets needed large 
ad hoc reduction factors for the type I migration rate as found in the isothermal approximation. 
In order to eliminate these factors, a simple semi- analytical model is presented that incorpo- 
rates recent results on the migration of low mass planets in non-isothermal disks. It allows for 
outward migration. The model is used to conduct planetary populations synthesis calculations. 
Two points with zero torque are found in the disks. Planets migrate both in- and outward to- 
wards these convergence zones. They could be important for accelerating planetary growth by 
concentrating matter in one point. We also find that the updated type I migration models allow 
the formation of both close-in low mass planets, but also of giant planets at large semimajor 
axes. The problem of too rapid migration is significantly mitigated. 
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1. Introduction 

The timescales of orbital migration of low mass planets as found in linear, isothermal 
type I migration models are very short for typical protoplanetary disk conditions. They 
are in particular shorter than typical growth timescales (Tanaka et al. 2002). This means 
that most protoplanets would fall into the star before reaching a mass allowing giant 
planet formation (~ 10 Mq). It is therefore not surprising that previous planetary pop- 
ulation synthesis models have found that isothermal type I rates need to be reduced by 
large factors (10-1000) in order to reproduce the observed semimajor axis distribution of 
extrasolar planets (Ida & Lin 2008, Mordasini et al. 2009b, Schlaufman et al. 2009). 

The necessity of such arbitrary reduction factors indicates a significant shortcoming in 
the understanding of the migration process. Additionally one finds that with universal 
reduction factors (independent of planet mass and distance), it is impossible to reproduce 
at the same time giant planets at several AU and close-in low mass planets (Mordasini 
et al. 2009b). This is inconsistent with observations (Howard et al. 2010). 

Several mechanisms were proposed that could slow down type I migration. An approach 
that recently gained significant attention was a more realistic description of the thermo- 
dynamics in the protoplanetary disk, in order to drop the simplification of isothermality 
(e.g. Paardekooper & Mellema 2006, Kley & Crida 2008, Kley et al. 2009, Paardekooper 
et al. 2010, Baruteau & Lin 2010). 

It was found that the migration rates (and even the direction of migration) in such 
more realistic models can be very diflFerent from the ones found in the isothermal limit. 
There are strong dependences on disk properties, like the temperature and the gas surface 
density gradient or the opacity, leading to diflFerent sub-regimes of type I migration. 
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Figure 1. Left panel: Schematic representation of the timescales relevant to the migration 
regime of a low mass planet. Trajectories of different gas parcels are shown in a frame of reference 
rotating with the planet (big black circle) around the star (center). Right panel: Specific torque 
as a function of planetary mass from the analytical model (lines, indicating different regimes), 
compared to the radiation-hydrodynamic simulations of Kley et al. (2009) (crosses). 

2. Updated type I migration model 

We present results of incorporating these improved type I migration rates into the 
planet formation code of Alibert et al. (2005). The new migration model, which will be 
presented in details in Dittkrist et al. (in prep.), includes the following mechanisms: 

Isothermal vs. adiabatic regime Gas parcels in the horseshoe region make a sharp U- 
turn close to the planet (Fig.[l] left). If during such a U-turn, the gas can cool quickly 
enough by radiation to equilibrate with the surrounding gas, the gas behaves in a locally 
isothermal way. On the other hand, if the U-turn timescale ru-turn is short compared 
to the cooling timescale Tcoob then the gas parcel keeps its entropy, and the process 
is approximately adiabatic. These two situations lead to different density distributions 
around the planet, which eventually translate into different torques. For both regimes we 
use the results of Paardekooper et al. (2010) for the migration rate. 

Saturated vs. unsaturated regime The horseshoe region only contains a finite reservoir of 
angular momentum, which can cause the torque originating from this region to disappear 
after a finite time. In order to check if such a saturation of the horseshoe drag occurs, we 
compare the libration timescale r^b and the viscous timescale 7~Yisc across the corotation 
region (Fig.[lj left). If Tvisc < 'T'lih sustained outward migration is possible. In the other 
case, Lindblad torques (plus some residual horseshoe drag) are acting, driving usually 
inward migration. In some cases, the (outward) migration of the planet itself can keep 
the corotation region auto-unsaturated in a feedback effect. 

Reduction of the gas surface density As the width of the horseshoe region increases 
with the planetary mass (e.g. Masset et al. 2006) and the viscous timescale increase with 
this width, saturation sets in when the planet reaches some mass, of order 10 M^. Such 
planets are particularly vulnerable to a rapid migration into the star, as the migration 
rate increases with planetary mass. On the other hand, starting gap formation reduces 
the gas surface density around the planet, reducing the torques. This is taken into account 
using the results of Crida & Morbidelh (2007). 

Gap opening criterion The transition to type II migration is built on the results of 
Crida et al. (2006). 

The right panel of Fig. [l] shows the specific torque (torque per mass unit) as a function 
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Figure 2. Left panel: Normalized torque as a function of semimajor axis at five times of the 
disk evolution. The vertical arrows indicate the locations of the convergence zones at 1 Myr. 
Right panel: The background shows the regions of inward (blue) and outward (green) migration 
in the disk. The lines show the migration tracks of seven planetary embryo inserted at different 
starting locations. Colors indicate migration regimes. Blue: locally isothermal. Green: adiabatic, 
unsaturated. Red: adiabatic, saturated. Brown: type II. 

of planetary mass. Positive torques correspond to outward migration. The line represents 
the result from our semi-analytical model, while crosses represent the 3D radiation- 
hydrodynamic simulations of Kley et al. (2009), indicating good agreement. 



3. Convergence zones and migration tracks 

With this model at hand, we have studied the migration of growing planets embedded 
in 1+lD standard alpha disk models (see Lyra et al. 2010). 

The left panel of Fig.|2] shows the normalized torque in the adiabatic, unsaturated 
regime as a function of semimajor axis and time for a typical disk. The torque depends 
on the gradients of the temperature and gas surface density which change because of 
opacity transitions. There are special zones where the torque vanishes, and where the 
torque gradient is negative. This means that inward of these points, migration is directed 
outward, and outward of them, it is directed inward, so that these locations are migration 
traps onto which migrating planets converge. 

The right panel shows the direction of migration as a function of time and semimajor 
axis. Green corresponds to outward, blue to inward migration. All protoplanets in the 
convergence zones migrate towards the stable zero torque locations. We find two conver- 
gence zones in agreement with Lyra et al. (2010), which seem to be a generic property. 
We note that the convergence zones are several AU wide, and therefore could concentrate 
a lot of matter in one point. This could have important implication for planetary growth, 
a process very recently studied by Sandor et al. (2011). We also see that the convergence 
zones themselves move inward on a viscous timescale (Paardekooper et al. 2010). 

The lines show migration tracks of seven embryos (initial mass 0.6 M^) migrating and 
growing in the disk. Each planet was simulated separately. The inner three protoplanets 
migrate outward to the inner convergence zone and remain attached to it for the rest of 
the evolution. They remain low mass planets (3 to 6 M^). The outer four protoplanets 
migrate both inward and outward to the outer convergence zone. They stay there until 
the corotation torque saturates because of the mass growth and fast inward migration 
sets in. Around 1 AU however, the planetary cores become so massive that gas runaway 
accretion is triggered, and the migration changes into the slow (planet dominated) type 
II regime. The final masses of these planets are between 4 to 6 Jupiter masses. 
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Figure 3. Concurrent growth and migration of protoplanets in the mass-distance plane. Plan- 
etary embryos start with an initial mass of 0.6 at different initial semimajor axes. Final 
positions are indicated with large black dots. Colors indicate different migration regimes as in 
Fig.|2] The left panel shows the nominal model, while for the simulation on the right, saturation 
is assumed to set in at a four times larger mass. 

4. Planetary population synthesis 

Figure |3] shows planetary formation tracks in the mass-distance plane, illustrating how 
planetary embryos concurrently grow and migrate using the undated type I migration 
model. Other settings and probability distributions are similar as in Mordasini et al. 
(2009a). Planets were arbitrarily stopped when they migrate to 0.1 AU. 

One sees that planets starting inside the inner convergence zone quickly migrate inward 
leading to close-in low mass planets. Further out, outward migration frequently occurs. 
Especially in the right panel, the imprint of the two convergence zones can be seen by 
two groups of giant planets. In the nominal model in the left panel, many embryos still 
migrate to 0.1 AU, but much less than with the full rate of Tanaka et al. (2002). It is 
clear that (giant) planet formation is no more suppressed by the loss of the embryos into 
the star. These results will be discussed in details in Dittkrist et al. (in prep.). 
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